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Abstract 

We consider asymptotic behavior of the correlation functions of the characteristic 
Oh, polynomials of the hermitian sample covariance matrices H n = n~ 1 A^ n n A m ^ ni where 

A m n is a m x n complex matrix with independent and identically distributed entries 
$ta a j and Qa a j. We show that for the correlation function of any even order the 
asymptotic behavior in the bulk and at the edge of the spectrum coincides with 
those for the Gaussian Unitary Ensemble up to a factor, depending only on the 

■ fourth moment of the common probability law of entries 9ia Q j, $$a a j, i.e. the higher 
^ ■ moments do not contribute to the above limit. 

in 

1 Introduction 

in 

Characteristic polynomials of random matrices have been actively studied in the last years. 
The interest was initially stimulated by the similarity between the asymptotic behavior of 
the moments of characteristic polynomials of a random matrix from the Circular Unitary 
Ensemble and the moments of the Riemann ^-function along its critical line (see |13|). But 
with the emerging connections to the quantum chaos, integrable systems, combinatorics, 
representation theory and others, it has become apparent that the characteristic polyno- 
mials of random matrices are also of independent interest. This motivates the studies of 
the moments of characteristic polynomials for other random matrix ensembles (see e.g. 

nu, na, a, ®, m, m : m, 0, ma, m. 

In this paper we consider the hermitian sample covariance ensembles with symmetric 
entries distributions, i.e. n x n random matrices of the form 

H n = n A* m ^ n A m ^ ni (1-1) 

where A m ^ n is an m x n complex matrix with independent and identically distributed 
entries ^ta a j and Qa a j such that 

E{a aj } = E{(a aj ) 2 } = 0, E{|a ai | 2 } = 1, a = 1, .., m, j = 1, .., n, , , 
E{^ 2l+1 a aj } = E{% 2l+1 a aj } = 0, I G N. 1 ' 
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We assume that m belongs to a sequence {m n \ c ^ =l such that 



m 



- -» c > 1, n — >■ oo. (1.3) 

n 



m,n— >co 



We below denote this limit as " lim 

Let A? 
(NCM) as 



Let Ai , . . . , An be the eigenvalues of H n . Define their Normalized Counting Measure 



N n (A) = %{\ { ; ) eA,j = l,..,n}/n, N n (R) = 1, (1.4) 

where A is an arbitrary interval of the real axis. The behavior of N n , as n — > oo, is studied 
well enough. In particular, it was shown in |17| that N n converges weakly in probability to 
a non-random measure N which is called the limiting NCM of the ensemble. The measure 
iV is absolutely continuous and its density p is given by the well-known Marchenko-Pastur 
law: 

p(A) = { -^>/(A + -A)(A-A_), A E a, ^ 
I 0, A £ a, 

where 

A ± = (l±v^) 2 , a = ((l-V~c) 2 ,(l + V~c) 2 ). (1.6) 
The mixed moments (or the correlation functions) of characteristic polynomials are 

F 2k (A) = JUdetiXj-H^P^dHn), (1.7) 



where is the space of positive definite hermitian n x n matrices, P n (dH n ) is a proba- 
bility law of the nxn random matrix H n , and A = {A-,-}^ are real or complex parameters 
that may depend on n. 

We are interested in the asymptotic behavior of ( j 1.7ft for matrices ( II. lft as m, n — > oo 
and for j = 1, ..,2k 

Ao + 6/np(A ), A G a, 



Aj_i A + ^/(n 7± ) 2 / 3 , A = A ± , {L * } 



where A± and o are defined in (11.61) . 

p is defined in ( 11 .5j) . and £ = {^j} 2 ^ are real parameters varying in [— M, M] C KL 

In the case of hermitian matrix models the asymptotic behavior of ( 1 1.7p was obtained 
by using the method of orthogonal polynomials (see [3j [19]). Unfortunately, the method 
of orthogonal polynomials can not be applied to the general case of the hermitian sample 
covariance ensembles (11.1ft - ( 11.2ft . In the paper [21] the method based on the Grassmann 
integration was developed to study the asymptotic behavior of the correlation functions of 
any any even number of the characteristic polynomials of the hermitian Wigner ensemble. 
Here we apply this method to the hermitian sample covariance ensembles (jl.ip - (jl.2p . 
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In [12] Kosters use the exponential generating function to study the second moment, 
i.e. the case k = 1 in (|1.7|) . It was shown that for Ao G o 

^-yF 2 (Ao + Ci/MAo)), Ao + 6/(np(A ))) = 2vrAr m C+ 1 / 2 

x e— m exp{nA + a(A„)(£i + 6) + 2k 4 } ^^^ (1 + o(l)), 



where 

A - c + 1 



A G a, 



a(A ) = <( 2A p(A ) " k 4 = - 3/4, (1.10) 

(1 ± v^)-^ 273 , A = A ± , 
7± is defined in (11.91) . and /14 is the fourth moment of the probability law of QWjfc, $tWjk- 
In P2] for the case c > 1, m = cn + o(n 1 / 3 ), k = 1 the asymptotic behavior at the edge 
of the spectrum (i.e. for Aq = A±) was also obtained: 



\t^ 2 (Ao + 6/(^7±) 2/3 , Ao + 6/(n 7± ) 2/3 ) = 2tt(1 ± ^cf 



(n— m) m+1/2 

x e 2 ^exp{n 1 / 3 a(A ± )(ei + £2) + 2/^(6, £ 2 )(1 + o(l)) 

with 

j4(x , !/)= Ai'WAi( j/ )-Ai M Ai-M 
z - y 

where Ai(x) is the Airy function 

Ai(x) = — / e is3/3+isx ds, (1.12) 

27r y 

5 = G C| argz = 7r/6 or arg z = 5ti/G}. 

In this paper we consider the general case k > 1 of (jl.7p for the random matrices ( II. ip . 
Define 

n(nW x . { MAo))- 1 i 71 2 (Ao + e/MAo)),A + e/(MAo))), A G a, 
Zy (c, An) — \ / \ (1.13) 

\ (n7 ± )- 2/3 F 2 (A + e/(^7 ± ) 2/3 ,A + e/(n7 ± ) 2/3 ), A = A ± , 

and denote 

2k 



J D 2fc (Ao) = nV jD(n) te,A ). (1.14) 
1=1 

The main results of the paper are the following two theorems: 

Theorem 1. Let the entries ^sa a j, 3la a j of the matrices U.l\) have a symmetric probability 
distribution with finite first 4k moments. Then we have for k > 1 



lim f„ n f\ \\l 2 n (\ \ F<2k ( A ° + ?/( n P(Ao)) 



cHk -i)/2 exp{k{k _ lK(c _ Ao + 1) V 1 }^ r s in(7r(&-& +J -)) , ( j , , , 



...,£ m )A(£ fe+1 , ...,£ 2 fc) I 7r(&-6 
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where F 2 k and p(X) are defined in U.7\) and M.5\) . A = (Ao, . . . , Ao) G M. 2k , X G a, 
£ = {£j}!=i, and and a are defined in $1. 10\) and $1.61) . 



Theorem 2. Let the entries ^sa a j, 9fta Q j of the matrices U.l\) have a symmetric probability 
distribution with finite first 4k moments, A = A± and let m belong to a sequence {m„}™ =1 
such that 

m n = cn + n 1//3 £ n , c > 1, (1-16) 
where e n — > 0, n — > 00. Tnen we nave /or > 1 



Sz to^5^o F » ( A ° + f/( " 7±)2/3 



n^oo (n7 ± ) 2fe / 3 _D 2 /c(A±) 

c fc(fc-i)/ 2 exp { 4 ^ _ !) K4 } r , & 



where F 2k and 7 ± are defined in U.7\) and M.9\) . A = (A±, . . . , A±) G IR 2A: ; £ = {£ 3 -}|=ij 
and K4 and A± are defined in U.10\) and U.6\) . 

The theorems show that the above limits for the mixed moments of the characteristic 
polynomials for random matrices (II. ip coincide with those for the Gaussian Unitary En- 
semble up to a factor depending only on the fourth moment of the common probability 
law of the entries a a j, i.e., that the higher moments of the law do not contribute to the 
above limit. This is a manifestation of the universality, that can be compared with the 
universality of the local bulk regime for Wigner matrices (see [7] and references therein). 

The paper is organized as follows. In Section 2 we obtain a convenient asymptotic 
integral representation for F 2 k, using the integration over the Grassmann variables and 
the Harish Chandra/Itzykson-Zuber formula for integrals over the unitary group. The 
method is similar to that of [21]. In Section 3 and 4 we prove Theorem [T] and [2j applying 
the steepest descent method to the integral representation. 

We denote by C, Ci, etc. various n- independent quantities below, which can be differ- 
ent in different formulas. 



2 The integral representation. 

In this section we obtain the integral representation for the mixed moments F 2 k (jl-7|) of 
the characteristic polynomials, i.e. we prove the following proposition 

Proposition 1. Let A 2k = A + £/(an) a , where A = diag{A , .., A }, £ = diag{£i, ..,£ 2 fc}, 
and 

p(A ), A G a, > 2 1 - 

7±, A = A±, 

" - < 2/3, x Cl, (2 - 2) 
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where a and \± are defined in U.6]) , and let F 2k (A 2k ) of ( II. 7]) be the correlation function 
of the characteristic polynomials. Then we have for every k 



D 2k \\ )F 2k (A 2k ) 



n 2fc: V-V) fe ( 2fc - 1 > 
2 k ir k e 2kn D 2k (\ ) 



2k 

IT 

3=1 



E 



v d &- 



2fe 



2fc ("1 _ , 

I! 1 w+2fc ; exp {2c m , n K 4 S 2 ((I - V)A ) f[ - 



vi 



i=i v i 

where V = diag{i>i, . . .,v 2k }, 



i=i 



vi 



1 + 0(n m,n —¥ oo, 



D 2k (\o) is defined in pTTS) and u is any closed contour encircling 0. 



(2.3) 



(2.4) 



To this end we use the integration over the Grassmann variables. The integration was 
introduced by Berezin and widely used in the physics literature (see e.g. [2] and [6J). For 
the reader convenience we give a brief outline of the techniques here. 



2.1 Grassmann integration 

Let us consider two sets of formal variables {/0j}?=D {V^j^u which satisfy the anticom- 
mutation conditions 

ipjipk + ipkipj = ^j*Pk + ^k^j = ipjipk + ^k^j = 0, Ji k = l,..,n. (2.5) 

These two sets of variables {■?/>.,' }™ =1 and {?/>,,• }™ = i generate the Grassmann algebra 21. 
Taking into account that if)? = 0, we have that all elements of 21 are polynomials of 
{V'jljLi and {^}" =1 . We can also define functions of the Grassmann variables. Let \ 
be an element of 21. For any analytical function / we mean by f(x) the element of 21 
obtained by substituting \ in the Taylor series of /. Since % is a polynomial of {%}" =1 , 
{^jYj=\i there exists such I that \ l = an d hence the series terminates after a finite 
number of terms and so f(x) £ 21- 

Note also that if \ is the sum of the products of even numbers of the Grassmann 
variables, then, according to the definition of the functions of the Grassmann variables, 
expanding (z — x) _1 into the series we obtain for any analytic function / 

/W " Z - /M, (2.6) 



n 



z — x 2?™ 



where Q is any closed contour encircling 0. 

Following Berezin we define the operation of integration with respect to the anti- 
commuting variables in a formal way: 

dipj = / dipj = 0, / ipjdipj = / ipjdipj = 1. (2.7) 



This definition can be extended on the general element of 21 by the linearity. A multiple 
integral is defined to be a repeated integral. The "differentials" dipj and dip k anticommute 
with each other and with the variables ipj and ip k . 
Thus, if 

k 

f(Xi, ■ ■ ■ , Xk) = ao + ^2 a hXh + a hhXhXh + ■■■ + a 1:2 ,...,kXi ■ ■ ■ Xk, 

Jl = l jl<32 

then 

f(Xi, ■ ■ ■ , Xk)d Xk---dxi = di,2,...,k- 



Let A be an ordinary hermitian matrix. The following Gaussian integral is well-known 



j,k=l j=l 



One of the important formulas of the Grassmann variables theory is the analog of this 
formula for the Grassmann algebra (see |2]): 

/n n 
exp | A 3\k^^k } IJ d $i = det A ( 2 - 9 ) 
A h. — 1 A — 1 



j,k=l j=l 



Besides, we have 

i 



/4 it it 

n^'p^ ex P { Yl A ii' c j''i-}Yl (l c j (I ''j = detAj 1 ,..,i, ; «i,..,*„ ( 2 - 10 ) 

p=l j,k=l j=l 

where A\ x \ ;si) .., s is a (n — q) x (n — q) minor of the matrix A without lines / 1; .., l q and 
columns s\, .., s g . 

2.2 Asymptotic integral representation for F2 

In this subsection we obtain (I2.3P for k = 1 by using the Grassmann integrals. This 
formula was obtained in [12] by using another method. We give here a detailed proof to 
show the basic ingredients of our techniques that will be elaborated in the next subsection 
to obtain the asymptotic integral representation of ( 11. 7ft for k > 1. 
Using ( 12. 9ft . we obtain from ( II. 7ft 

2 n _ 

F 2 (A 2 ) = El / e-i»«=i d* 2>n 

^ 2/n \ / n \ 2 n 

= E [] e ^ P=1 9 " e3=1 P " <**2,J (2.11) 



r /• ™ * / i™_ _ \ea s e? ps 1 

1 / 1111 I 1 ~ ~ ^P^qWprWqr )e s=1 p=1 dw 2| j, 



E 



ct=l r=l p,q=l 
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since for any a = 1, .., m and any r = 1, 2 we have according to ( 12. 5 j) 

(n v 2 / n \ 2 

E^V) =f5Z a «9^) =°- ( 2 - 12 ) 
p=l ' ^ g=l ' 

Here {"0ij}"^=i are t ne Grassmann variables (n variables for each determinant in (jl.7p ) 
and 

s I 

d^ s ,i = Y[Y[d^ jr d^ jr . (2.13) 

r=l j=l 

In view of 01 .2p and (12. 5 p we get 

f 2 / 1 n — \ ] i 2 ^ " _ 

E< ]^[ ( 1 E U <xpa aq l/J pr l/Jqr ) f = 1 E E E {^°P a "9} V>VV (2.14) 

r=l ^ p,q=l ' * r=l p,q=l 

^ n n 

"^^2 ^ ] ^ ] ^ {_ a ap 1 a aq 1 a ap2 a aq2} ' l Pp 1 l' l Pq±l' l Pp22 l Pq2'2 
pi ,i}i=l P2,g2=l 

l 2 ™— 1 — — 1__ 

= 1 ~ ~ E E ^Mfcr + ^ E ^plMtf^tf ~ ^ E ^pl^M^l 
r=l p=l p^q p^q 

+ — E ^pi^v^ P 2^p2 = det Q 2 n) + — ^ ^ 1ppllppl1pp2' l Pp2, 



p=i p=i 



where and are the matrix with Grassmann entries 

*P = {E^*y > gP = l-n- 1 ^, (2.15) 
I p= i J r,t=l 

/i 4 is the 4-th moment of the probability law of ^sa a j, ^ta a j of (11. 2p . and k 4 is defined in 
Thus, (l2~iTj) and (12TT4D yield 



^ 2 (A 2 ) = / e tr ^' A2 (detgW + ^ E^pi^)^^ 
E Q^ 1 / e tr ^ A2 det--Qr ) (E^i^ p2 ^ 2 ) 9 ^ 2 



^ E Q (^^ / ef ^'drf-gf^d * 2 ,_ 



(2.16) 



>A /m\ n! (2/t 4 ) g 
: ^ ™ 2g ^ 



Here we used the symmetry of ijj lp , ipi p and 



4 , P i^pi4 , p2^p2f(^pi^pi^p2^p2)d^ p id4> p idilj p2 dilj p 2 = /(0, 0,0,0). 
To compute ^.g we use the following lemma 



Lemma 1. Let A be any p x p matrix and let I be a positive integer. Then we have 

r Ar AU 

det l A = K p , l I —g-dp(U), (2.17) 



where 

p-i 



^=(-i) p(p - i)/2 5 P i ri( z + s ) ! ' s p=U si ( 2 - 18 ) 



s=0 



U is a unitary matrix with eigenvalues {uj}? =1 , W is a matrix which diagonalizes U and 

dn{U) = A 2 ( Ml) . . . , u p )d wf[^ (2.19) 

where duj means the integration over the circle u = {z : \z\ = 1}, dW is the Haar 
measure over the unitary group U{p), and A(u\, . . . , u p ) is the Vandermonde determinant 

of Uj-S. 

Remark 1. 1. Lemma[l\ is a particular case of the superbosonization formula which was 
proved in the physics paper [16]. We give below (see Subsection 2.4) a different proof for 
this simple case. 

2. Since both sides of $2.17\) are analytic functions of a^j, we can take A with not 
necessary complex but also with even Grassmann elements. 

3. Combining A2.17\) and 1(2.9) we get that for any p x p matrix A 

r a ti au 

tr 



e* AVV 'd9rt = K P>1 I -—^-dfMiU), (2.20) 



det p+( U 

where iff® = {^ s =i ipsjipsr}^ r= i an d d^ff Pi i is defined in \2.13) 
Using Lemma [T] and Remark 1.3, we obtain 

„tr A 2 *("- 9) +trw2-n- 1 tr w 2 ^ q) 



h q = K 2 „ J det m - q+2 W 2 d* 2 , n - q dti(W 2 ) 

r g trw 2 +tr A 2 J7 2 -n- 1 tr U2W2 
= V,V, J ^^^^ KUMW.) (2.21) 

= K 2 , n - q J d»(U 2 ), 

where U 2 and W 2 are unitary 2x2 matrices, and dfi(U 2 ), d fi(W 2 ) are defined in (I2.19p . 

Recall that we are interested in A2 = Ao,2 + £2/ (na) 13 , where Ao,2 = diag{Ao,Ao}, 
fa = diag{£i,£ 2 }, and a, are defined in <|27T]> . (Q . Substituting (12191) in ([2721]) and 
using that functions det(J — n -1 ^), tr A U 2 , and det U 2 are unitary invariant, we obtain 
from ( 127211) 

h, q = K 2 , n _ q [ I e tr Ao,^ + (na)-^tr hwvw -Q ( X r_|)T! {V1 _ V2 fdv(W)^l2.22) 
J J r=1 v r {Am) 

U) 

= K l2 n zl f f tr A , 2 y 2 +M-^tr bwnw TT i 1 ~ v r) m ~ q , _ y d (w] dvidv 2 

n 2(n- q) Jf e 11 v n- q+ 2 ^ ^^^(2^)2' 



where u is any closed contour encircling 0. The integral over the unitary group U (2) can be 
computed using the Harish Chandra/Itsykson-Zuber formula (see e.g. [H], Appendix 5): 

Proposition 2. Let A be the normal p x p matrix with distinct eigenvalues {aj}f =1 and 
B = diag{&i, . . . , b p }. Then for any symmetric function f(B) of {bj}^ =1 we have 

J f e^ AU * BU A 2 (B)f(B)dUdB = S p J *Wf(b 1 ,...,b p )dB, (2.23) 

U( P ) 

p 

where S p is defined in Ii2.18\) . dB = Y\dbj, dU is the Haar measure over the uni- 

3=1 

tary group U(n) and A(A), A(S) are the Vandermonde determinants of the eigenvalues 
{ai}U>{h}UofAandB. 

This and formula ( I2.22p yields 

_ 2n^ l a^K 2 ^ q f tr Aq aVa+n i-e a -pti £ 2 v 2 A (1 - v r ) m ' q (vi - v 2 )dv 1 dv 2 

2 ' g ~ n^-i) f e l\ v r q+2 (ei-6)(2«) 2 - 



Hence, since 



"' {n - q + mn -^ = 2nn^(l + 0(l/n)) 



(n — q)\ n 2n q 

we get (ESD for k = 1 from (|2TT6|) . fl2~T8|) . and ( I2T24D . 

2.3 Asymptotic integral representation for i 7 ^ 

Using (EU) and fl2TT2|) . we obtain from flTT]) (cf. (j2TTT| ) 

2fc n _ 
f /" E E ( X l- H )p,q^ prVV 

F 2k (A 2k ) = E\ / er=i*.«=i d* 2fc 



E W eS=1 P " II' P=1 q " d ^2k,n\ (2.25) 



2fc n 



a=l 
m 2fc 



{n T — * i T — * i i Tib Tl 

f L^L fps^ps -i-r -i-r / 1 \ - _ — \ 



a=l r=l P,9=l 

In view of (II. 2p similarly to (I2.14p we get 
E< JJ (l a ap a aq i) pr ip gr ^ I = detQ^f 



P>9= 

+ ^ E det(Q£ ) )^'^ sl ' S2 ^^ 1 ^ 1 ^ 2 V' PS2 +n- 2 $(*), (2.26) 

(l<i2,Sl<S2 P=l 
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where Q 2 ™ } is defined in (1215)1 . det(Q$?) ( ' 1,ia;ai,aa) is ( 2fc - 2) X (2Jfe - 2) minor of matrix 
Qzk without lines si, s 2 and columns Zi, Z 2 , «4 is defined in (ll.lOp and is a polynomial 
of the variables {(^ _1 ^2fe )r-,«}r s=i an d 

^ n q 

n ' 1(X j^ = -^YI^pi^psj, l = (h,---,l q ), s = (sx,...,s q ). (2.27) 
p=i 3=1 

Now we use 



Lemma 2. Set A = {aij}^-!, 6 = ^ ere l,s is defined in $2.27\) . Let <& r (A,b) be 

an analytic function of the variables and {bj^} and let (1 — e)n < r < n, < I < en 

with some sufficiently small e > 0. Then there exist an absolute constants Cq,C% such 
that 

Mn-^ln-'a^li^d^^ < C max \$ r (A,b)\ ■ [ Ji% (*)d* 2fc , r , 

^(*) = e^ )A »det"-'gW. (2.28) 



where 



The proof of Lemma [2] is given in Subsection 2.4. 

Denote the expression multiplied by k 4 in the r.h.s. of (I2.26P by n~ 1 X. Write 



(det Q% + + n- 2 $ (¥) 

m ! 

fci!fc 2 !(m- fci - fc 2 )! V~~^ 2fc y Vn 



It is easy to see that the terms in (12. 29ft such that k\ + k 2 > em give the contribution of 
order e~ enlogn and thus can be omitted. Hence, (j2.25p . (I2.26p . and Lemma [2] yield 

F 2k (A 2k ) = (1 + Oin- 1 )) J d* 2k , n e tr *- A - 

/ " n \ m (2 30) 

(detQ^ + ^1 £ det(Q(?)^^--)^^ l ^ 1 ^ 2 ^ 2 ) , 

where and are defined in (I2.15p . 

To compute the r.h.s. of (I2.30p we use the Newton binomial formula and observe 

that the term with p s = p, in the product ]J q j=1 (// 1 V,,. /; , '.^/^ ) can 

be expressed in terms of (I2.27P with an additional factor n~ l . Therefore, according to 
Lemma [2] it suffices to consider only the terms with p s ^ pi or, taking into account the 
symmetry, the term p\ = n, p 2 = n — 1, . . . , p q = n — q + 1 with coefficient n\/(n — q)\. 
Thus, we can write 

MM = (1 + 0(»-)) ± (™) {n H :J 2K J' h k „, (2.31) 
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where 



hk, q = J (2.32) 

p $<*) = II ( E det(Q£ ) ) a? »^ sf ' sf) ^^K^f^f)' ( 2 - 33 ) 
where J^ q is defined in (12.28j) . Note that 

det(Q^ ) ) (h, ' 2;Sl,S2) = det(Q^" <?) ) (il ' i2;Sl ' S2) + 

are polynomials of {^is}^^) {V^'sl^s^i with the sum of the coefficients of order 
O(q), m, n — > oo. Hence, using Lemma [2j we get 

hk, q = J fltl^P^W^kAl + 0(q/n)) =: J 2fcj9 (l + 0(q/n)), (2.34) 

According to Lemma Q] and ( 12. lOf) we can rewrite l2k, q as 



e tTK 2k ^ k - q) +tlQ^- q) V 



det m-g+2k y 



n 2k 



2 k 



X 



X 



/ II Wd^pid^Pi^ 1 P — Q+1 (2-35) 

/• e trA 2fc *<7 9) + tr Q2r 9) ^ 

^2fc,m-g y dfi(V)d#2fc,»-9 ^-1+^ y 

( det(g(r 9) ) ( ' 1 ' 2;Sl ' S2) det(A 2fe -n-V) ( , 1 , 2;si , S2) ] 9 . 



• h<h,si<S2 

Besides, the Cauchy-Binet formula (see [9]) yields for 2k x 2k matrices A, B 



detA^ s ^B {luh , sl ^ = \f- 2 det(x - AB) 



2dx 2 

h<h,si<S2 



x=0 
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Thus, using again Lemma [T] and Remark 1.3, we obtain 

r e trA 2fc <- 9) +trQ<r 9 V 
hk, q = K 2k ,m- q J dfi(V)d^2k,n- q detm _ g+2fc v 

FT AlL. [ FT K 2k,ldKW s ) tTQ^- q) (A 2k -n-W)±W p -±z p tTW p 

l\:2mz*J 11 det 2k+l W s 6 [ ' ' 

s=l * s=l " 

r p tTA 2k U+tT(l-n~ 1 U)V 

K 2 k,m-qK 2 k,n-q / d fi(V)d fi(U)- 
J I 



det m- g +2fc y det n-q+2k jj 

dz s J ftr 2feil d//(W s ) e tr(i-n- 1 c/)(A 2fc -n- 1 y) EW p -^Z p tTW p 



n j n det 2fc+i ^ 

s=l s s=l 



e tr A 2fc c/ det — «(! _ det m -"(l -n-^W, 



i 

! 

= ^ fc ,n- 9 y d et n -* +2fc U ~ 

11 2^ 11 det 2fc+1 ^ 

s=l * s=l A 



Besides, 



9 , 



det^l-rT^Wp) = e ' * =1 "(1 + 0(g/n)). 
p=i 

Substituting this to (I2.36P and using (12 . 1 7p . we get 

e^^^det^^l-n- 1 ^) 



hk, q = (1 + 0{q/n))K 2 k,n- q J 



det n-q+2k u 
1 I 



1 2tt^ 3 1 A det 2fc+i w 

s=l s s=l A 



[1 + 0{q/n))K 2Kn . q I detW _ g+ L ^|((1 - r i - 1 C/)A 2 ,)^(C/) 



Recall that we are interested in A2& = Ao,2fc + £,/{ na Y \ where Ao = diag{Ao, .., Ao}, 
£ = diagj^i, ..,^2fe}) an d a 5 /3 are defined in 02. ip . (I2.2p . Thus, 



hk, q = (1 + O(n^) + 0{q/n))K 2k , n _ q , +2fc 



e trA 2fc t/ det m-^ 1 _ n -i [/ ) 

det n- g +2fc ^ 

x S|((l - n- 1 U)A )dfi(U), (2.38) 



Let us change variables to U — W*VW, where If is a unitary 2k x 2k matrix and 
V = diagj^i, . . . , v 2 k}- Since functions det(J — n' 1 U), S 2 ((I — n~ 1 U)A ), and det U are 
unitary invariant, (12.381) implies 

2k 

hk, q = (1 + 0{n-f>) + 0(q/n))K 2k ^ q I f[ pL f d ^W)e trw * VWA ^A 2 (V) 

I 3=1 J (2.39) 

det — _ n -iy) 

det ,- g+2fcy ^ f((I-n ^)Ao). 
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where u is any closed contour encircling 0. The integral over the unitary group U (2k) can 
be computed using the Harish Chandra/Itsykson-Zuber formula (I2.23p . Shifting Vi — > nvi, 
we obtain 



l2k,q 



n 



2k(n-q) ]J 



2k 



n 

3=1 



2ni 



2k (\ _ v \m-q 

n 1 n - L S!(I -V)(l + 0(n-P) + 0(q/n)). 



Hence, since 



1=1 v 



n—q+2k 
I 



ris=o 1 ( ri -^+ s ) 



(27T)V fc2 e- 2fcn (l + 0(l/n)), 



(n-q)\ n 2k(n-q)+q 

we get 023} from ( I2T341 . and 



(2.40) 



2.4 Proofs of Lemmas [D, [2] 

Proof of Lemma [T] Let A be a normal matrix. Then we can set A = VqAqVq and 
U = W*UoW, where A = diag (ai, .., a p ), Uo = diag (ui, .-,u p ) and Vo, W are the matrices 
diagonalizing A and [/ correspondingly. We obtain 



/ := 



Atau 



det p+l U 



dU 



e tr W ™ A 2 (lil) _ )ttj) ) 



i=i 



p+« 



dii(W)J[ 

3=1 



duj 



Shifting integration with respect to W as WVq — > W and using ( 12. 23ft . we obtain 



p 

ajUj 

i k ei=1 A ( M i' • • • ' u p> TT 

,=q >$ ^^n 2 „ 

3=1 



, M p ) T-r d Uj 



Qp 



A(A C 



det 



e ajUj duj 



a(a ) n «s 

3=1 

1 P>P-1 



p+l - s 2%i 



A(A ] 



det 



J j,s=l,0 



p+l-s-l "I P'P" 1 



(p + i-s-1)! 



P n(p + Z-s-l)!A(A)) 

s=0 



i,s=i,o 

[-l)^\detU 

u(p+i-s-i)\ ii( P +i- s -i)\ 

s=0 s=0 



and (I2.17P is proved for the normal A. 

Let now A be an arbitrary matrix. According to the polar decomposition, we can write 
A = SW, where W is a unitary p x p matrix and S is a diagonal p x p matrix. Since we 
proved ( 12. 17ft for any normal A, we proved it for S = diag{e* Ql , . . . , e iap }, a\, .., a p G JR.. 
Besides, it is easy to see that both sides of (12. 17ft is analytic functions of the elements of 
S. Therefore, we are proved (I2.17P for any A. 
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□ 

Proof of Lemma [2] According to Lemma [1] and (12. 28 p . we have 

f f f gtrAafctfW+trtl-n-^WjV 
/ ^ 2 l i md^ 2k)r = K 2k ^ l J J det m-l+2ky dfl(V)d^ 2k>r 

f e t " v det r (A 2k -n- 1 V) ^ (2Al) 

= K ^, m -l J det m-l+2ky d ^ V ) =■ 1 

It is proved below (see Section 3 and 4 (note that if we change v — > Aq(1 — v) in J we 
obtain the integral like in (13. 3p and (14. ip )) that 



2fe 2fe 



i j i>^Srf ef=i ^ n^-^nA^nr-^w, (2.42) 

~ 3=1 3=1 1 J ' 

where a and (5 are defined in (12. ip and ( 12.2D and 

£ = L G C : |*| = (^iA V /2 }. (2.43) 



n 

Moreover the integral outside of the any n-independent neighborhood of v± = (Ao + — 

1) /2±7rAop(Ao) give contribution 0(e~ nC ), hence we can can deform co near z = f^^Ao) 
such that \v — Ao| > 5 on u. Thus, if we define 

((•••)> = r 1 J(...)g il (*)d* a , r (2.44) 

the definition is correct. 
Using ( 12. 6p . we get 



(I a, , T-r <l Ih 



($ f .) := ($ P (n-i*w )ff «)> =f* r {A,b) n n 



2vri 



/ 2fc 1 TT 1 

x ( . n : „_i,,T,(rh n r - 1 , i / • 



(2.45) 



Thus, to prove Lemma [21 we have to estimate the expectation above. Expanding the 
functions into the series with respect to {^ 2k )i,j}, {<Jj^}, we get 

II , II 



2k 2k 



ee e (n(»-(»s^n(»-w n «d"-w 0* 

ij,7,s'i,j=l%,a =1 *J =1 J,s 1 >3=1 Z,s 

:=EEE tf (fcMyn»riic r ' 
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To estimate the moments {M({l it j}, {tj ■§})}, we introduce the generating function 

F(C,z) := (expjn-HrC^ + n- 1 ^^}), (2.47) 

l,s 

where ( = {Cij}ij=i- It is easy to see that the derivatives F((,z) with respect to {Cij} 
and will give us the moments {M({l it j}, {tj ■§})}. 

Using Lemma [Hand then integrating over d^2k,r, we obtain 

tr V+tr (A^+n-^-n^V^+n- 1 £ zj-ai^ 

nc,z) = ^f^ J - det ^ v — — wmmw 

k r «try 



J 



where J is defined in (I2.4ip . Moreover, according to (12. 9 j) - (I2.10|) . $i(A 2 fc— n -1 V, rir 1 ^, n~ l z) 
is a polynomial of the entries of A 2 fc — rT x V and of {Ci,j}j {^is} with n-independent coeffi- 
cients and degree at most 2k and such that the degree of each variable in $i(V, ^ _1 C> n~ 1 z) 
is at most one. Here we also used that the integral over d^2k,r can be factorized in 
{Vv^m}8=i- Besides ; 

^(n-V.n^Cn-^) = det(A 2fe - rrV) + /(A 2fc - n-^n-^n" 1 ^ (2.48) 

where /(A 2 fc — n~ 1 V, n' 1 ^, n~ l z) contains all terms of $i which includes {Ci.j} or { z j-§\- 
Recall that we are interested in A = Ao + ^/(na) 13 , where Ao = diag{Ao, .., Ao}, £ = 
diag{£i, . . . , £ 2m }, and a and j3 are defined in (12. ip and (12. 2p . Change the variables 
Vj — > nvj, j = 1, ..,2k, where {vj} are the eigenvalues of V, and replace the integration 
over the unit circle by the integration over u of (12. 43ft 

jy- f „ntr V 

F ^ Z ) = j . gg^j j detm - l+ 2 kv ^ k -V,n~\,n~h)d^V) (2.49) 
We have from the description of $x and (I2.48P 

(2.50) 

where C(V) > is bounded for G u with a; of ( 12.43ft (recall that we can deform u near 

1 /2 

z = (^^Ao) such that |i> — Ao| > 5 on ui ). Since {M({Zjj}, {tj 3 })} are the derivatives 
of the generating function, we can write 



IfilAa - V.n-'C.n- 1 *)! < C|det(A - V)| f[ (l + ^Cul) JT (l + ^^Lil), 



, ■ f .'p . J ^1 y 



This, (12T49|) . and (12~50|) yield 



2fc 



|M({/ M },{^})| < min ^-!e c l*l-^^Wn rin fyle ®-****® (2.52) 

ij=l 
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Choose Q id = {C G C : |C| = k,j/C}, S T _ = {z e C : \z\ = tj-jC}. Then (1X521) yields 



2 /v 



\ M ({kj}i{\s})\ < n v^c^nvHi * 



2A- 



(2.53) 



M =1 J,5 



Thus, if loijl > C and |&^| > C in (1X46]) we obtain Lemma M from (I2H]L (g35D - (12T46D . 

□ 



3 Asymptotic analysis in the bulk of the spectrum. 

In this section we prove Theorem (U passing to the limit m, n — > oo in (I2.40P for Aj = 
\ + €j/np(\ ) : where p is defined in (TL5|) . A G ct with a of (FO|) . and ^ G [— M, M] C R, 
j = l,..,2k. 

To this end consider the function 

V(v, A ) = -A f - c m , n log(l - u) + logw + S 1 *, (3.1) 

where 

^ ri* ^0 C m n -|- 1 C m n C m n 11 . » 

C -=n' * = 2 + ^ l0g ^-2 l0g V (3 ' 2) 

Then (JXJ) and (J22QD yield 

W»K . . . , w 2fe ) J] tZv^l + o(l)), (3.3) 

.7 = 1 

where L^fc is defined in (11. 14ft . 



UIQ 



A (0 i=i W i (3.4) 

2fc 



x exp {2c m , n /£ 4 5 2 ((/ - 7)Ao) II i=^}> 



and 

2fe 

,2 . -2fcre 4 -a(Ao) £ Cj 

n fe p(A ) fe(fe 1] e ' =1 
Z2fc ~ 2 2 % 2fe c fc / 2 ' ^ ^ 

Now we need the following lemma 

Lemma 3. The function dtV(v, Xq) for v = A 1 ^ 2 e lv , <p G (— 7T, 7r] attains its minimum at 
v = v ± := A 1/2 e^ := A °-^" + 1 ± mp(A ). (3.6) 
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Moreover, if ip ^ U n (±ipo) := (±y9 — w 1 / 2 logn, ±y? + n 1 / 2 logn), £/ien we have for 
sufficiently big n 

^(A 1/2 e^A )>^^. (3.7) 

n 

Proof. Note that for ip G (— 7T, 7r] 

KV(A 1/2 e^, A ) = -Aj /2 cos ^ - ^ log (l + A 1 - 2A 1/2 cos y?) + log Aq 1 + S*, (3.8) 
where S* and c mi „ are defined in (13. 2p . Thus 

^ TOT// \ -1/2 w \ \ \l/2 • C m n /A 

—&V{\ Q e^,A =V sinyjfl , 

d V v l + A 1 -2A 1/2 cos^- 

^XV^e* A ) = Af cos <p (l i C ^( A _° 1/2 ) (3.9) 

d V 2 V 1 + Aq 1 - 2A 1/2 cos^ / 



2c mn sin y?/A 



o 



(1 + Aq 1 -2A 1/2 cos^) 2 ' 

-I Icy 

and </9 = ±<y2o of (13. 6ft are the minimum points of HXV{\q ' e^, A ). Writing 

\-l/2 g j n 

V± := V(-u±, A ) = =RA 1/2 sin ip ± «Vo ± iCm,n arcsin —2 — | , (3.10) 

1 + A — 2A cos <^o 

we conclude that 

3?V>±,A ) = 0. 

Expanding KV(Aq 1/2 e^, A ) into the Tayl or series and using (I3.9P - (I3.10p . we obtain for 

if G U n [±(p ): 

^(Ao 1/2 e^,A ) = ( 7r V(Ao)) 2 ( ^ Tyo)2 + 0( ^3/ 2lQg 3 n); (3 n) 

where <p is defined in (13.61) . This and c m ^ n — > c, m, n — )■ oo imply for ip G" U n (±tpo) 

^(Ao V2 e- A)>^^. 

n 

The lemma is proved. □ 

— 1/2 

Note that \Vj\ — A , j — l,..,2k. Since £i, ...,£2*: are distinct, the inequality 
| A (T) /A (0| < d and (E2D yield 



/n n 2 Al 

® .. <b W n (vi, . . . ,v 2 k) Y^dvj 

where 



< Cm k e 



u = {z G C : \z\ = A 1/2 }, (3.12) 
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W n and Z 2 k are defined in (13. 4p and (13. 5p respectively, and 

U± = {ipe (— 7T, 7r] : \±(p -<p\ <n" 1/2 logn}, (3.13) 
U v ,± = {z = Xo 1/2 e l ^eU ± } 

with ip of (13. 6p . 

Note that we have for ip e C7± in view of ( 13. ip and ( I3.10P as m, n — > oo 

V(A - 1/2 e-, A ) = V ± + (± - (i^j) tfe^&^p?- + f±((p =F ?o), (3.14) 

where /±(<p =F y?o) = =F V^o) 3 )- Shifting ^ =p V^o — >• f° r e ^± an d using ( 13 . 1 0[) 

we obtain 



^(Ao) „ rA , „ , H 2k-i)/2^sr f °''^>^ Ag) 

(l/n)°'" 



S 71C_L_ 2Al ?7.C q 2fc 

n e -^^-(2fc-i)^ J+W ) fi e -T^^-(»-i)i(w-w>) ndiPj=rrT l)a , 
1=1 j=h-i i=i s=1 a , 

(3.15) 

where a s = {ctj} 2 ^ is a permutation of s pluses and 2k — s minuses, 

d±((pj) = —= , V s = diag{e i(vi+vo) , .., e i(Vs+w) , e^ 1 "^, ^ e i(^-vo)^ ( 316 ) 

vAop(Ao) 

G a (<p u ..., V2k ) = 2c m , w « 4 g 2 ((Z - V')Ao) J] \ _ 1/2 - ~ : [ 1 , 

f=i 1 - A o ' e'^+w) 1 - A 1,2 e<vr-vo) 



2k-s 

II 

1=1 r=s+l 



c± = Qj- ~ " m ; n ±)2 ) A - 1 e ±2 ^, C/ n = (-n- 1 ^ log n,n-^ log n), (3.17) 

Define 

s 



f ^EUw) n 
I s = e ]_]_(<^ - tpi)dv a [tp u . . . ,ip s ) 

det (e^ 9fe Vi _1 V ■ • • , P.) (3-18) 



OO 



pi,..,p s =0, 



- - /; 
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where di> s ((f\, . . . , tp s ) is a measure on Q n )S := ( — logn, \ogn) s which is symmetric in 
((fi,..,(p s ) and g(^p) is a function such that g(ip) = Cip(l + o(l)), n — > oo. Note that 
if we take the term of (I3.18P such that p Sl = p S2 , Si 7^ then this term is zero since 
d ZAf((/?i, . . . , ip s ) is symmetric in (ipi, . . . , ip s ). Moreover, the order of 

det {(n-^gi^^/p^l^ 

is n~( pi+ " +Ps ^ 2 and if {pi, ..,p s } 7^ {0, 1, .., s — 1} the order is less than rTT 5 * 5-1 )/ 2 . Hence, 
denoting by the sum over all permutations {pi, of {0, 1, .., s — 1}, we obtain 

/ nfe^)) PiA ^ 1 '"'^) £i ^^ 1 '---'^)( 1 + ( 1 )) 

-n.7- </ „_i 



S 



,rf iti^-i •, / n^) i " iA (^-'^) d ^(^'---'^)( i +°( 1 )) ( 3 - 19 ) 

n s{s-l)/2 [\ j=0 j\ J 



A(c/((^i), .., g(ip a ))A(<Pi, ••, f,(y>i, ...,<pa)(l + o(l)), 



'n.s 



Since g((p) = C(p(l + o(l)), n — >■ 00, we get 

I= nS ( s -i)/m""',! / A 2 (^,..,^)d^(^,...,^)(l + o(l)). (3.20) 

o 

Consider T a of ( 13. 15f) with ai = . . . = a s = +, a s+ \ = . . . = a 2 k = — . Since the function 

2c min K 4 5 , 2((/ - Vk)A ) EI 1 _°-i/2 ei(yi+c ,^ o) is symmetric in (y>i,..,p a ) and ((p s+1 , ..,(p 2k ), 

changing variables as y/n(fj — > cpj and using formulas ( 13.181) - f l3.2Qf) , and formula for the 
Selberg integral (see, e.g., [14], Chapter 17), we obtain 



lo Z n 2k 



T 

- 1 a 



llj=oJ-llz=o / 1=1 7=1 



c +¥ > 2 



2k 



l=s+l 



xA% s+1 ,..,^) II e"^(l + o(l)) = -^§^-^(1 + 0(1)), (3.21) 



where Co )S (£) is an n-independent constant. This expression is of order 0(1) for s = k, 
and it is of order o(l) for s 7^ /c. Hence, only the terms of (I3.15P with exactly k of {a.,} 2 ^ 
pluses contribute in the limit (jl,15p . If we take s = k we obtain 

exp{z7r(^i + ... + £ fc - £ fc+i - ... - 6fc)} c fc(fc-i) K4 (c-A +i) 2 e - 1 ^ 3 22 \ 

rij=i rL=fc+i(£i _ 6) 

A fc2 (2mp(A )) fc2 cfc(fc _ 1)K4(c _ Ao+1)2e - 1 e-^+-+^-^ +1 --^) 
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Hence, since it is easy to check that 

4vr 2 A2p(A ) 



c+c_ 

we get from ( I3.2ip and ( 13 .22 j) that T a of (I3.15P with a\ = ... = ctk = +, otk+i = ... = 
a2k = — has the form 

Mk+i) p i7r(6+-+a-5fc+i---6ft) N9 . 

(2zvr)^ n (6 - 



In view of the identity 



df i sin(7r(^-6 + 0) | fc det /e-fe-^) - e *tt*+i-fc) ^ ' 



A(£i, .., a)A(6+i, », 6fc) (2i7r) fc A(ei, .., &)A(&+i, ..,&*) 

2k 

the determinant in the l.h.s. of f 1 3 . 2 3 j) is a linear combination of exp{iir ^2 otj^j} over 

3=1 

the collection {aj}^ x , in which m elements are pluses, and the rest are minuses. By the 
virtue of the following formula (see j 18 j . Problem 7.3) 

2 ~^Tk — ; 7\ =det{(o i -6i) } . (3.24) 



n,-,i=i(%-6. 



the coefficient of exp{i7r(£ fc+ i + ... + £ 2 fc - 6 _ ••• _ 6fc)} is 

det {(&+,- fe)" 1 }^, (-1)^ 



Other coefficients can be computed analogously. Thus, restricting the sum in (13 . 1 5[) to 
that over the collection {Qtj}|*j_, in which exactly k elements are pluses, and k are minuses, 
and using fl3.23[) . we obtain Theorem [1] after a certain algebra. 

4 Asymptotic analysis at the edge of the spectrum. 

Let now A = A + (for Ao = A_ the proof is similar) and Xj = \ + + ^j/(n'j + ) 2 ^ 3 , j = 1, .., 2k, 
where A + and 7+ are defined in ( 11. Qh and (11. 5ft . and £1, . . . , C,2k £ [ — M, M] C M. 
According to ( 12. 3 j) we have 

1 2k 

( ^ +) ( 2 A fc t/ ) 3 i? 2 fc (A 2fc ) = Z 2k jw n {v u v 2k ) Y[d Vj (l + o(l)), (4.1) 
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where D 2 k is defined in fU.14p . 

W n (vi, . . .,v 2k ) = exp | - n^V {+) (^) + ^ 2/3 ~ Vl + n ( c m,« ~~ c ) JZ 10 ^ 1 ~ 



X 



7 

z=i i=i 1 i=i 



1 } ^ {2c m , n . 4 ^ 2 ((/ - v)Ao) n n ^ 



(4.2) 



A(0 I ' t=i l - Vl) U"i 

V {+ \v) = -X v - clog(l - u) +bgu - S+, (4.3) 
5+ = -1 - - clog(l - (1 + v^T 1 ) - log(l + v^, (4.4) 



and 



_ -n 1 / 3 a(A + ) J] £ 3 +2fc(nc-m)log(l-A~ 1/2 ) 

Z 2fc = L 2fc e > =1 , (4.5) 

n k(2k+l)/3^2k(k~l)/3 e -2kK, 4 -2k(nc-m) log(l-A^ 1/2 ) 

L 2fc = 2 2 % 2fc c fc / 2 ' ^ 4 ' 6 ^ 

We need the following lemma 

Lemma 4. T7ie function Ry' + '(t)) /or u = A+^e 2 ^, y> G (— 7r,7r] attains its minimum at 

v := A; 1/2 = (1 + Vc)- 1 . (4.7) 

Moreover, if v G Wo = {f G C : f = A + 1 ^ 2 e* v , </3 G (— 7T, 7r]} ; |t> — t> | > 5, where 5 is small 
enough, then we have for sufficiently big n 

$V {+ \v) > C5\ (4.8) 

Proof. Similarly to (13.81) - (13. 9p we have 

^(+)( Uo ) = __s R y(+)( Uo ) = __5Ry(+)( Uo ) = 0, (4.9) 

(/ $tV {+ \v ) = 6. (4.10) 



dip 



Hence, ip = is a minimum point of the function ( A + e 1 *' 3 ) , and ;JH^ + )(A_ ) _ 1 ^ 2 e l¥ ') 

is monotone increasing function for ip G [0, it) and monotone decreasing function for 

If G (-7T,0]. 

Expanding 3?W +) (A; 1/2 e^) into the Taylor series we obtain for \ip\ < 5 similarly to 

icmi: 

3?V (+) (A; 1/2 e^) = ^ 4 /4 + 0(p 5 ). (4.11) 
This and monotonicity of (A + 1 ^ 2 e %Lp ) for 99 7^ imply for 9? G (— 7r, 7r], |^| > 5 

Since \v - v | = 2A; 1/2 | sin(^/2)| < A; 1/2 |y?|, we get ( HTTP . □ 
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Note that \vj\ = A + 1 ^ 2 , j = l,..,2k and according to ( I1.16P c mi „ — c = o(n 2 / 3 ), 
m,n — > oo. Since £i,...,£2& are distinct, the inequality |A(T)/A(£)| < G\ and (14. 8 \ 
yield 



7 = 1 

UJq\Us(vo) ^0 ^0 J 



Z 2fc J j) .. j> W n (v 1 ,...,v 2 k) Y[dvj 

< Cl n k ^ k+1 ^ 3 e- c ^ 1+0 ^ +Csnl/3 , m,n -> oo, 



where 

t^Oo) = e w : |u - vol < S}. 

Hence, 

^ k ^] 3 F 2k (A 2k ) = Z 2k j W n (v 1 ,...,v 2k )f[dv J (l + o(l)) + 0(e- Cn ). (4.12) 
Since 

^+)(« b ) = ^y(« b ) = o, 

we have for | v — Vq \ < 5 

y (+)(„) = 7 - 2 ( w - Wo ) 3 /3 + 0((« - t; ) 4 ), |u - v \ -> 0. (4.13) 

Thus, we can write for v satisfying |u — u | < 5 

V< + )(«)= 7 -V(«)A (4.14) 

where x(t> ) is analytic in the 5-neighborhood of i>o with the analytic inverse z((p) (we 
choose x( v ) such that x(t>) el for wet). 

Changing variables to = z(tpj), j = 1, .., 2/c, we rewrite (I4.12p as 

, 2k 2k 1/3, 2fc i / \ 



(n7 + ) 2fc2 / 3 



2 Cm , nK4 s 2 ((yr-z)A ) n i^jA(2) ^\ zfUfi) , ,„ /1%N r „ s , A „s 

1=1 ^n^) rf ^ i +°( i ))+ o ( e ) ( 4 - is ) 

2fc 



=:L 2fc / % b .., fe )[[^ 3 (l + o(l)) + 0(e 

~ i=i 

where L 2k is defined in (14. 5ft . 

Z = diag{z(<pi),...,z(<p 2 fc)}> (4.16) 
^ = {^C: z(^) G £/f(v )}. (4.17) 
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Moreover, we have from ( I4.14p 

X(v ) = 0, -f- z X(vo) = 1, (4.18) 

hence 

< d < \ X '(v)\ < C 2 , \v - v \ < 5. (4.19) 

If a = {z G C : \z — Zq u \ < 5}, then x(^) i s a closed curve encircling (p = and 
lying between the circles o"i = {ip G C : |<^| = C\5} and cr 2 = W G C : |y?| = 6*2(5} for 
< d < C 2 . We have from (HT8D 

z(0) = «b, 2/(0) = 1, < Ci < \z'{ip)\ <C 2 , if G (4.20) 

According to Lemma IU W/^ + \v) > for t> G Us(vq) and we get Ky?^ > for y>j G Z/j^, 
i.e., 

cos(3argv? i ) > 0, (pj e U SiV , 
where Ug iV> is defined in (I4.17p . Hence, Ug^ can be located only in the sectors 

— 7r/6 < axgip < 7r/6, 7r/2 < argy? < 57r/6, 7tt/6 < aigip < 3n/2. 

Besides, x is conformal in a (see (I4.19P ). hence angle-preserving. Taking into account 
that xi v ) 6 K for d G M, the angle between Uq and the real axis at the point v is 7r/2, 
and that t/^ v is a continuous curve, we obtain that U$ }ip can be located only in the sectors 

tt/2 < arg^ < 5tt/6, 7tt/6 < argy? < 3tt/2. (4.21) 

Note that we can take any curve U((p) instead of Ug i9 provided that U(cp) and ujq \ Us(vo) 
are "glued", i.e., the union of z(U((p)) and u>q \ U$(vq) form a closed contour encircling 0. 
Let us take 

U(<p) = {if G C : argip = 2tt/3, y? G 

U {y? G C : argy? = 4vr/3, y? G U f/ lj5 U U 2 , s , 

where a = {v G C : |i> — Vo\ < 5}, is a curve along x(^) from the point of intersection 
of the ray argy> = 27r/3 and x{9&) to the point ipij of intersection of £7$^, and x{9&) 
(7r/2 < argy?!^ < 5?r/6), and f/ 2 ,<5 is a curve along x{@&) from the point of intersection 
of the ray argy> = 47r/3 and x(9a) to the point <p 2) s of intersection of Us iV> and x{9&) 
(7it/Q < argy? 2 ,5 < 3n/2). According to LemmaHand ( j4.14j) , 5 = r 3 cos3y3 > C > 0, 
where r = |y?i i( s|, y?o = argy9i i( 5. Since < C\ < r < C 2 , we have 

cos 3(p Q > C/Cl > 0. 

Moreover, it is easy to see that cos(3 arg ipi) > cos 3y2o along Ui t g (since cos 3x is monotone 
increasing for x G [7r/2,27r/3] and monotone decreasing for x G [2tt/3, 57r/6]). This and 
\<fj\ > C\ imply for cpj G L^g 

8 > C > 0, </?i G f/ 1>5 . (4.22) 
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Also we have from (I4.20p 

\z(<Pj) - v \ < C 2 \<Pj\ < C, (pj e x(p)- 
This, fl422|) . c m ,„ - c = o(n~ 2 / 3 ), m,n -> oo, and (TOO]) yield 



< e 



-Cn+o(n) 



V?i e U h5 , (pj e Ufa), j > 1. (4.23) 



Hence, the integral over t/i^ does not contribute to the l.h.s. of (14. 15f) The same statement 
we can prove for L^s- Thus, we have shown that integral over Us }ip in (I4.15P can be replaced 
to the integral over the contour 

T= {(p G C : arg<^ = 2tt/3, <p G X (*)} U {(/? G C : arg <p = 4vr/3, up G (4-24) 

According to the choice of Z, we have 

%J = rJ, ^gT, (4.25) 

where r 3 - = |<^-|. 
Set now 

a n = {<p G C : M < e,- 1 ^- 1 / 3 }, 

where e n is defined in (II . 1 tip . Note that we assume that w -1 / 3 — > 0, n — > oo. In other 
case we can take o n = {uo G C : \(p\ < logn/n 1 / 3 } and the proof will be similarly. 

It is easy to see that a n C x(&) f° r sufficiently big n. Besides, we have from f )4.20p for 

<f> G a n 

z{cp) = v + uj + O^n- 2 / 3 ), n -»■ oo, 
^) = 1 + 0(en 1/2 n- 1/3 ), n ^ oo. 1 ' 

Taking into account (051) . (OB)) . fl4T25|) . and 



, 1 - 

log 


< 


w - *(<p) 


1 - u 




1 - U 



wc 



obtain for <pi G Z \ cr n , ^ G /, j = 2, .., 2/c 



W (if x,..., <f 2k) 



(4.27) 



where ri = > Sn^^n 1 / 3 . Since ra 1 / 3 r 1 > en. 1 ^ 2 for ^ £l \ cr n , the integral over Z \ er r , 

-3/2 



is 0(e Cen ) as n — > oo. Hence 



(n7 + ) 2fc2 / 3 



2 A: 



F 2fc (A 2fc ) = L 2fc (l + o(l)) / W(vj u ...,f 2 k) Y[duj J + 0(e' C£ " y2 ), (4.28) 



where L 2k and W^i, • • ■ , f2k) are defined in flJE} and (Q5l) . This, (E3D, and ( )4~26l) 
imply 

e 



T 2fc(2fc-l)« 4 

2fc p a ^ 2fc 



(n7 + ) 2fc2 / 3 



*-2fcJ 



,4fc 2 



-n 7 " 2 £ ^t|^ ; A($) 



1=1 1=1 T 



2k 



A(0 



(4.29) 



x(l + <%i, . . . ,v?2fc)) EI dl Pj + °( e 



-Ce, 



-3/2 
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where 5(ipi, . . . , <f2k) collects the reminder terms which appear when we replace z(tpj) — > 
v + <Pj + O(^), z'M 1 + 0{ Vj ) and log ^ + 0(cp§), j = 1, .., 2k. Hence 

\6(<fi, . . .,(f2k)\ < C(\(pi\ + . . . + \(f2k\)- 

Changing variables in (I4.29P as r y~ 2 / 3 n 1 / 3 Lfj — >■ icpj we obtain in new variables 

. . . , ip 2k )\ = \8[i^n- l ' Z ip u . . . ,i 7 2/3 n- 1/ V2fc)| < Cn~ l l\\^\ + ... + |^|). 
Therefore, using (II. 9p . (14. 5p . and (14. 7p . we obtain 

D^{K) „ fA , L 2k e^ 2k -^(l + o(n-V 3 )) 



(n7+) 2fc2 / 3 V ' " V 4fc 2 (_i 7 -2/3 n l/3)fe(2fc+l) 

2fc 2k 2k 



x / e ; =i ;=1 — d<pi + 0(e Ce " ) 



x 



e 4fc(A ; -l) K4( _ 1) ^ c fc(fc-l)/2 (1 + ( n -V3)) (4.30) 

, 2fc 2fc ok 



where S is defined in ( 11. 12j> . 
Consider 



2k 2k nu 



/ det l^e^V^m I (4.31) 



(27T) 2fc 

Integrating by parts, we have for j > 3 



j,i=i j=i 



dipt 

\{j - 3)^'- 4 + ^rv v?/3+ ^^- 



Applying this identity to each line, starting from third, we observe that the first term in 
the r.h.s. gives zero contribution. Repeating this procedure and replacing and rearranging 
the lines, we obtain from (I4.3ip 



■k r ( -\ 2k 2k 

K(0 = J^J (-l) fc(fc - 1)/2 det |^epe^ /3+l ^} t U d ^ ( 4 - 32 ) 
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where j = qjk + rj, qj = 0, 1, rj = 0, 1, .., k — 1. Thus, 

fc 

Kg) i" r v Afe„ ...QJAfe,, ..,( ik )Un. t g^g^,» 

A; 

C_l-)fc(fc+l)/2 TT „, 2k 2k 

e i e ^ !=i n<%> 



(27T) 



2/v 



5 n<-<i* ]"J ]"J _ ^) sign (t _ Js ) i=l 

_ (4.33) 

where the sum is over all collections 1 < ji < ... < Jfc < and A , . . , £ Jfc ) is the 
Vandermonde determinant of with j ^ j 1; .., Consider 



U- 2 



Oet < > 2k 2k 91, 

J A(6,..,e,)Afe + i,..,6,) L\ ^ r { } 



(2vr) 2fc J A(ei,..,6)A(a+i,..,6., 

fc 

According to the identity (13.241) the coefficient of f| in (14.341) is 

s=l 

^fc /_-|\fc(fc+l)/2 



(27r)2fc fc 

II II Oj s - ¥>t) sign (t - j s ) 
s=1 t^=ji,..,jk 

Thus, K(£)/A(£) is equal to (14.341) . and (I4.30j) yield the assertion of Theorem [2j 
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